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ABSTRACT 


This  report  deals  with  the  study  of  electrostatic  Troves  in  bouadeo 
hot  plasmas  and,  except  for  a  few  dc  measurements,  is  theoretical. 

The  study  is  broadly  divided  into  two  parts.  The  first  section  In¬ 
vestigates  .?aves  in  a  diode  system  with  no  dc  magnetic  field  applied.  A 
dc  analysis  of  the  equilibrium  plasma  produced  within  a  thermionic  diode 
is  undertaken  to  provide  a  basis  for  the  rf  analysis,  and  expected  plasma 
densities  in  the  diode  are  computed  for  a  variety  of  practical  situations. 
Electrostatic  wave  resonances  in  the  diode  are  predicted  by  using  the 
hydrodynamic  model  for  the  norunifcrm  plasma.  The  impedance  of  a  uniform 
plasma  diode  is  obtained  by  using  a  kinetic  model  of  the  plasna.  This 
model  enables  us  to  take  into  account  the  end  plate  electron  absorption 
loss,  a  process  similar  to  that  which  causes  end  plate  diffusion  in  the 
Q-machine.  The  absorption  loss  is  fcund  to  have  a  large  effect  on  the 
impedance  of  the  diode. 

The  second  part  of  the  report  deals  with  the  study  of  guided  waves 
along  a  cylindrical  column  of  Maxwellian  plasma  in  a  magnetic  field.  A 
dc  study  of  the  plasma  column  is  first  conducted.  Theoretical  density 
and  current  profiles  are  obtained  and  are  compared  with  the  nes  •.  '•#>  re¬ 
sults.  Since  an  rf  analysis  using  a  self-consistent  dc  solution  i?  '.e 
involved,  the  plasma  column  is  approximated  by  a  uniform  plasma  with  a 
sharp  boundary  and  no  drift.  To  obtain  the  rf  fields  in  such  a  column, 
a  plane-wave  solution  for  an  infinite  Maxwellian  plasma  in  an  applied 
magnetic  field  is  first  obtained.  By  superimposing  infinitely  many 
plane  waves,  the  fields  within  the  column  are  constructed,  and  matching 
the  appropriate  fields  at  the  boundary,  a  dispersion  relation  is  derived. 

The  solutions  to  the  dispersion  relation  reveal  the  existence  of  a 

new  type  of  unstable  waves.  When  only  the  electrons  are  asrumed  to  respond 

to  the  electric  field,  the  surface  waves  which  propagate  when  _  < 

ce  pe 

are  found  to  be  unstable.  When  the  ion  motion  is  included,  additional 

unstable  surface  waves  ire  obtained  whenever  _  ^  .  A  studv  of  the 

ci  pi 

instability  shows  that  it  is  due  to  finite  Laraor  radii  effects  and  is 
driven  by  the  transverse  energy  of  the  particles. 


ill 


While  the  "electron”  surface  *ive  instability  is  relatively  week 

and  can  easily  be  stabilised  by  nakine  _  >  „  .  the  ”ioa"  surface 

ce  pe 

wave  instability  is  fouod  to  be  very  stroi^,  and  the  stabilizing  condition 

of  *  ,  >  -  is  not  easily  attained  in  hieh  density  plasms  such  as 
ci  pi 

those  in  fusion  Machines .  Thus,  tbe  unstable  ion  surface  saves  asy  have 
a  serious  effect  on  tbe  containaeat  of  tbe  fusion  plazas. 
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I.  INTRODUCTION 


The  object  cf  this  study  is  an  investigation  of  electrostatic  waves 
in  bounded  hot  plssaas.  While  the  bulk  of  prior  work  on  electrostatic 
waves  is  restricted  to  infinitely  extended  plasmas,  our  study  emphasizes 
the  effects  due  to  finite  plasma  boundary.  Originally  we  hoped  to  avoid 
complications  due  to  Instabilities  by  restricting  our  work  to  near  equi¬ 
librium  plasmas.  Yet  in  the  study  of  a  perfectly  Maxwellian  magneto¬ 
plasma  column,  we  uiscovered  unstable  surface  waves  which  are  due  to 
finite  Larmor  radii  effects. 

Plasmas  in  perfect  thermostatic  equilibrium  are  generated  inside  a 
black  body,  and  our  work  on  the  subject  of  equilibrium  plasma  generation 
and  the  waves  in  such  plasmas  is  summarized  in  Section  A  below.  A  near 
equilibrium  column  of  plasma  is  obtained  in  a  magnetic  field  between 
electron  and  ion  emitting  end  plates.  Prior  work  on  electrostatic  waves 
in  such  plasma  columns  and  our  own  work  which  lead  to  the  discovery  of 
unstable  waves  is  summarized  in  Section  B.  A  brief  review  of  kncwn  in¬ 
stabilities  in  Section  C  places  the  newly  discovered  instability  into 
proper  perspective. 

A.  PLASMAS  IN  PERFECT  THERMOSTATIC  EQUILIBRIUM 

Theory  and  interpretation  of  experimental  data  on  plasmas  often 
suffer  because  the  deviations  from  equilibrium  are  not  known  exactly. 

A  situation  in  which  all  of  the  constituents  are  in  equilibrium  may  be 
found  inside  a  black  body.  A  study  of  the  dc  states  of  an  equilibrium 
plasma  bounded  by  thermionlcally  emitting  walls  is  undertaken  in  Chapter 
II.  An  earlier  work  by  Nottingham1  considers  the  diode  geometry  con¬ 
sisting  of  two  parallel  infinite  walls  emitting  only  electrons.  In  such 
a  system  the  electron  density  within  the  diode  spuce  was  found  to  be 
highly  nonuniform  and  became  very  low  at  a  distance  of  a  few  Debye  lengths 
away  from  the  wall.  We  show  that  under  equilibrium  conditions  a  small 
amount  of  ions  emitted  from  a  thermionic  emitter  can  neutralize  the 
electron  space  charge  to  form  a  plasma  in  the  diode  space.  This  is  il¬ 
lustrated  by  computations  based  on  diode  walls  made  from  pure  refractory 
metals.  The  effects  of  enhanced  ion  emission  from  the  walls  due  to  the 


1 


introduction  of  alkali  vapor  in  the  diode  space  are  also  examined.  Some 
experimental  measurements  for  a  sodium  plasma  are  given  and  compared  with 
the  theoretical  prediction. 

Electrostatic  waves  within  the  equilibrium  plasma  diode  are  studied 
in  Chapter  III.  The  nonuniform  sheath  effects  are  investigated  by  means 
of  the  hydrodynamic  model  of  the  plasma.  Previously  the  hydrodynamic 

2  ?  4 

model  was  u^ed  by  Weissglas  and  later  by  Parker^  to  explain  the  Dattner 
resonances  in  a  nonuniform  plasma  column.  As  opposed  to  these  previous 
works  on  free  boundary  sheaths,  our  study  is  directed  to  the  cathode 
sheaths  bounded  by  an  emitter. 

In  the  second  part  of  Chapter  III  we  present  a  study  of  end  plate 
loss  effects  by  using  a  kinetic  model  of  the  plasma.  The  end  plate 
losses  are  due  to  the  same  mechanism  which  causes  end  plate  diffusion 
in  a  plasma  column  such  as  that  produced  in  the  Q-machine.  When  par¬ 
ticles  iti  the  magnetized  plasma  column  are  absorbed  by  the  end  plate 
and  re-emitted  with  different  guiding  centers,  the  net  result  is  a 
diffusion  of  the  particles  away  from  the  column.  In  the  diode  system 
consisting  of  infinite  parallel  planes,  we  are  not  concerned  with  dif¬ 
fusion,  but  the  process  of  absorption  and  re-emission  causes  rf  energy 
loss  in  the  diode.  An  electron  emitted  from  the  wall  acquires  rf  energy 

as  it  traverses  the  diode  space,  and  the  entire  energy  is  lost  as  it  is 

6 

absorbed  at  the  other  wall.  Previous  work  oy  Hall  used  kinetic  theory 
in  considering  a  plasma  capacitor  but  assumed  perfect  reflection  of  the 
particles  at  the  boundary  and  thereby  neglected  the  end  plate  loss  effects. 
In  our  study  we  include  the  absorption  effect  at  the  end  plate,  but 
because  of  the  extreme  complexity  of  kinetic  analysis,  we  were  able  to 
study  only  the  case  of  a  uniform  plasma  in  the  diode  space. 

B.  ELECTROSTATIC  WAVES  IN  A  PLASMA  COLUMN 

In  a  practical  plasma  the  emitting  plates  are  not  infinite  in  extent. 

In  fact,  plasmas  are  often  produced  between  two  circular  emitting  plates 
which  are  separated  by  a  distance  several  times  greater  than  the  plate 
radius.  Confinement  of  the  plasma  in  such  a  column  is  achieved  by  ap¬ 
plying  a  magnetic  field  along  the  axis.  A  dc  study  of  such  a  column  of 
magnetized  plasma  is  given  in  Chapter  IV.  Theoretical  density  and 


2 


current  profiles  are  computed  from  an  appropriate  solution  to  the  Boltz¬ 
mann  equation  and  a  comparison  is  made  with  some  experimental  results. 
The  study  indicates  that  the  plasma  column  may  be  approximated  by  a 
uniform  plasma  with  a  sharp  boundary.  To  simplify  the  study  of  waves  in 

the  plasma  column,  it  is  essential  to  assume  that  the  column  is  uniform. 

7 

In  I957  Gould  and  Trivelpiece  used  the  cold  plasma  model  to  predict 
electrostatic  wave  propagation  along  a  uniform  plasma  column  such  as  the 
one  shown  in  Fig.  1.1.  The  dispersion  relations  for  cylindrically  sym- 


metic  waves  along  such  a  column  are  shown  in  Fig.  1.2.  When  03  >  a>  , 

c  p 

propagation  is  confined  within  two  frequency  ranges,  0  <  a>  <  03  and 

H 2  2  P 

0)c  <  a)  <  +  u) Whether  the  plasma  completely  fills  the  waveguide 

or  not,  these  passbands  are  unaltered  as  long  as  a>  >  03  .  However, 

c  p 

when  a)  <  03  ,  filled  and  unfilled  waveguides  have  different  passbands. 
c  p 

For  the  completely  filled  waveguide,  the  two  passbands  are  0  <  03  <  co_ 

r 2  2  c 

and  03  <  03  <  /03  +  03  .  When  the  plasma  does  not  fill  the  waveguide, 

P  v  P  c  /  /  2  2  \ 

the  upper  passband  is  as  before  1 03  <  03  <  /  03  +03  f,  but  the  lowest 


but  the  lowest 


order  mode  of  the  low-frequency  band  propagates  in  the  band  0  <  03  < 

7  2  p 

(o3p  +  Ojp/2  as  shown  by  the  dotted  line.  This  particular  mode  is  a 


FIG.  1.1.  CONFIGURATION  OF  PLASMA  IN  A 
CYLINDRICAL  WAVEGUIDE. 
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FIG.  1.2.  DISPERSION  RELATION  FOR  WAVES  IN  CYLINDRICAL  PLASMA 
COLUMN.  Dispersion  relations  for  waves  in  a  cylindrical 
column  of  cold  plasma  are  given  for  two  cases.  The  diagram 
on  the  right  shows  dispersion  relations  with  u)c  =  1.2  fflp 
and  on  the  left  dispersion  relations  with  o>c  =  0.5  <»p  are 
shown.  The  quantity  ka  is  the  product  of  wavenumber  and 
the  waveguide  radius. 

/  2  2 

surface  wave  in  the  frequency  range  <  a>  <  J  (cd^  +  coc)/2  in  that  the 
fields  are  strongest  near  the  edge  of  the  plasma.  Such  a  surface  mode 
cannot  be  excited  in  a  filled  waveguide.  In  Chapter  VI  attention  is 
directed  toward  this  surface  wave  which  becomes  unstable  when  the 

temperature  effects  are  included  properly. 

8 

Agdur  and  Weissglas  used  the  hydrodynamic  model  of  the  plasma  to 

o 

point  out  the  temperature  effects.  Later  Jayson  and  Lichtenberg'  used 
an  electron  plasma  model  with  Maxwellian  velocity  distribution  along  the 
magnetic  field  but  with  zero  transverse  temperature  to  compute  dispersion 
relations  for  waves  in  filled  waveguides.  Kuehl  et  al1^  used  a  fully 
thermal  plasma  model  but  with  infinite  magnetic  field  to  compute  the  dis¬ 
persion  relation  in  the  identical  geometry.  The  models  of  Lichtenberg 
and  Jayson  and  Kuehl  et  al  yield  Landau  damping  of  the  modes,  while 

tne  hydrodynamic  model  does  not.  The  study  of  waves  in  fully  thermal 

11  12 

unbounded  magnetoplasma  has  been  conducted  by  several  authors,  ’  but 
the  waves  in  a  magnetoplasma  column  have  not  been  treated  before  due  to 
the  extreme  computational  complexity  involved.  This  paper  reports  the 
study  of  waves  in  a  magnetoplasma  column  and  the  discovery  of  surface 
wave  instabilities. 
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In  preparation  for  the  study  of  eloctrostatic-wave  propagation  in  a 
aapttlztd  plasma  column,  a  dielectric  tensor  appropriate  for  an  infi¬ 
nitely  extended  Maxwellian  plaswa  is  obtained  in  Chapter  V.  Our  dielec- 
trie  tensor  is  stapler  than  the  "dielectric"  tensor  presented  by  Stix.  J 
The  simplicity  of  our  dielectric  tensor  results  from  subdividing  the 
plasm*  current  into  polarization  and  magnetization  currents.  A  cylin¬ 
drical  wav*  solution  a  ipropriate  for  a  uniform  column  of  Maxwellian 
slasaa  is  synthesized  from  plane  waves  and  the  dispersion  relation  for 
the  waves  is  obtained  in  Chapter  VI.  At  low  magnetic  fields  (o>c  <  o^) 
the  surface  wsve  along  the  column  of  a  Maxwellian  electron  plasma  is 

shown  to  be  unstable,  the  nature  of  the  instability  is  studied  by  means 

14  15 

of  Derfler’s  stability  criteria.  *  The  surface  waves  which  occur 

near  the  ior.  cvclotron  frequency  when  oi  .  <  to  .  are  also  found  to  be 

ci  pi 

unstable.  Since  the  complete  analysis  is  so  vast,  only  the  lowest-order 
cylindrically  symmetric  mode  is  numerically  analyzed. 

C.  INSTABILITIES 

Plasma  instabilities  can  be  divided  into  two  general  groups,  micro¬ 
instability  and  mscroinstsbility.  The  instabilities  which  cannot  be 
derived  from  the  standard  magnetohydro dynamic  equations  but  which  depend 
on  detailed  microscopic  equations  for  the  plasma  are  classified  as  micro- 
instabilities.  We  are  interested  primarily  in  the  microinstabilities 
which  occur  in  a  finite  geometry. 

Among  the  macroinstabilities  in  a  current  carrying  plasma  are  the 
well-known  kink  and  sausage  instabilities  «’hich  are  caused  by  the  inter¬ 
action  of  the  plasma  with  its  seif-magnetic  field.  Another  instability 

which  occurs  in  a  plasma  with  a  directed  current  is  the  Kadomtsev  insta¬ 
ll  17 

bility.  The  Kelvin-Helmholtz  instability  occurs  in  a  stratified 

fluid  in  which  the  adjacent  layers  are  in  relative  motion. 

The  above  instabilities  can  be  obtained  from  the  fluid  model  of  the 
plasms.  The  second  class  of  instabilities,  the  microinstabilities,  may 
be  divided  into  two  broad  groups.  The  first  are  those  instabilities 
which  are  caused  by  the  peculiarities  of  the  velocity  distribution  and 
may  be  called  velocity  space  instabilities.  The  second  group  of  insta¬ 
bilities  depends  on  the  plasma  geometry.  A  simple  example  of  the  first 
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ii.  dc  suns  or  wouxm  canon  plash* 


Crftee  ^ltnu  pr&Jacrt  it  tfc«  Ubontcrj  differ  f roe  •  true  tberno- 
<y»mc»I  equilibria*  plaana.  Ia  suck  noaequillbrlua  plasmas  tlte  energy 
CitnkaUaa  of  tkc  ooettitatau  (electros,  ioa>,  neutrals  and  radiation) 
deals  tee  significantly  tram  lolUasM'i  probability  lee  (exp  -  E/KT)  . 

These  <a*l«Ucss  ere  an  apparently  mark  sot  title  source  for  the  discovery 
of  sea  tnstseiil tis* ,  datifsef  by  nature  to  uitmately  restore  equillbriun. 
Theory  and  lataryrtituoe  of  experimental  data  frequently  suffer  because 
the  donations  Iron  equilibria*  are  aerer  kaoen  exsctly.  A  situation 
a^ara  all  ccontitasats  are  as  close  as  possible  to  strict  themodynaadc 
squill  brass  any  be  found  si  this  a  boundary  consisting  of  vails  which 
mars:  sea  sally  «!  electrons,  ions,  and  neutrals .  le  this  chapter  me 
*1.12  aun  the  dc  properties  of  nqsullbrlaa  plasans  bounded  by  such 
therussraceiiy  actnsf  sails.  Waxen  cal  results  are  girea  for  sone 
practical  situations. 


a.  Tusu  non 

lx  this  sect: ox  se  isrostigste  theoretically  tae  plasma  produced 
unseen  the  t so  initials  plane  nutting  walls  as  shone  in  rig.  c.  1.  There 
is  ac  net  dc  poteailiei  across  lbs  diode,  and  the  emitter*  eatii  electrons 
and  ions  ix  none  last  metier  ratio  ~£  3  *{».;•*  (L)f  which  Is  fixnd 
t.  tan  j  eS'isuc  considers!  loos  of  the  eaten  si. 
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In  the  past  a  great  deal  of  work  has  been  done  on  diodes  consisting 
of  thermionic  emitters.  The  majority  of  these  papers  deal  with  the  study 
of  dc  and  oscillating  states  with  a  nonzero  dc  potential  applied  across 
the  diode.  An  early  work  on  the  thermionic  diode  with  zero  net  potential 
was  given  by  Nottingham,1  who  treated  the  problem  in  the  absence  of  ions. 
Without  the  ions  he  found  that  the  electron  density  was  highly  nonuniform, 
ana  ii  the  spacing  was  much  greater  than  the  Debye  length,  the  density 
at  the  aldplsns  was  very  low.  Eichenbaum  and  Hernqvist2^  included  the 
ion  emission  but  studied  the  problem  using  the  collisionless  Boltzmann 
equation,  never  allowing  the  empty  regions  of  the  phase  space  to  *ill  up. 
Figure  2.2  illustrates  the  model  used  by  Eichenbaum  and  Hernqvist  with 
tje  shaded  region  of  the  ion  phase  space  unfilled.  Such  a  nonequilibrium 
analysis  is  not  realistic  when  there  is  no  net  potential  across  the  diode. 
Langmuir  examined  an  equilibrium  case  but  with  one  of  the  planes  at 
Infinity. 


IONS  ELECTRONS 


FIG.  2.2.  PHASE  SPACE  SKETCH  FOR  THE  ELECTRONS  AND  IONS. 

In  our  problem,  we  assume  a  strict  thermodynamic  equilibrium  and  fill 
the  phase  space  completely  according  to  Boltzmann's  probability  law,  and 
hence  the  scale-height  law  is  used  for  the  electron  and  ion  densities. 

By  taking  the  potential  to  be  zero  at  the  midplane  (x  =  0),  Poisson’s 
equation  becomes 

=  "  F*  ^n+(x)  “  n_(x)3 

=  -  r  {n+(°)  exp[-  7f]  -  n»  exp[S]} 


(2.1) 


At  this  point  normalized  quantities 


5 


> 


and 


r\ 


eV 

A 


a 


n+(o) 


with  the  Debye  length  at  x  =  0, 

SqKT  \l/2 

n  (0)  e2/ 


are  introduced,  and  Eq.  (2..1)  is  thus  reduced  to 


(2.2) 


Let  us  first  consider  the  case  of  a  <  1,  the  "electron"  rich  case. 
The  boundary  condition  at  §  =  0  is 


T1  =  0  ,  iy  =  C  . 

By  introducing  a  new  variable 

v?  =  e11  ,  (2.3) 

the  solution  to  Eq.  (2.2)  is  found  to  be 

F(a,v)  =  K(a)  -  §  (2.4) 

where  F(a,w)  is  the  elliptic  integral  of  the  first  kin d2^  with 
modulus  a, 
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Cl 
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•ad  1(a)  x  p(a,l )  is  the  couplete  tlllytlc  integral.  T*kl«  far  tfcwt 

2r 

Integrals  arc  available. 

By  Inverting  Bq.  (<:.*).  a*  obtain 


»  *  »»>  (a;  - 

•here  sa[E(a)  -  5A  i»  *be  Jacobi  elliptic  fsacUoa1  ntk  an  Sr  V  i 
a  and  argument  1(a)- i.  Beace, 


-  =  -2  -ia  sbj(2)  - 


»  -» 


Equation  (2.7/  enables  os  tc  cxlrala’e  '.be  petaactial  as  a  Iocum 
of  z  provided  a  is  kaon.  Generally,  a  IS  *  qazxxity  «•  seek  ** 
terns  of  a  sialiar  quantity  glees  at  the  eadtter  surface.  Let 


.  c 


le  knot  from  tbs  scale- height  laa  that 


n_(L)  s  s_(:)  expi-^, 

and 

n+(L)  «  n+f.)  expr--^  . 

Hence, 

n+(L>  V-i 

=  ^rnr  exp;~£v 

or 


zc  =  a' 


FIG.  2.3.  RATIO  OF  ELECTRON  AND  ION  DENSITY, 

Of2  =  n  (0)/n_(0) ,  IN  MID  -PLANE  AS  A  FUNCTION  OF 
SEPARATION  DISTANCE  2L  IN  UNITS  OF  DEBYE  LENGTH 
A  FOR  DIFFERENT  VALUES  OF  p2  =  n  (L)/n  (l)  . 


For  a  approaching  unity  we  have 

12  2 

sn(u,a)  «  tanh  u  +  -jj  (1-0!  )  sech  u[sinh  u  cosh  u  -  u]  ...  (2.12) 


and 


K(a)  w  la 


L(i-a) 


1/2J  *  * 


(2.13) 


For  small  values  of  w  or  yf$/OL,  we  also  have 

F(a,w)  m  w  -  -jr  w2(i+a^)  ... 


(2.x':) 


By  shifting  the  origin  to  the  wall. 


5'  -  -  5  +  L  , 


using  Eq.  (2.9)  and  letting  the  other  wall  go  to  infinity  (a-*  l),  we 
obtain  from  Eq.  (2.7) 

r\  «  -  2  -In  sn[F(l,/p)  +  1] 

»  -  2  -to  sn[  Jb  +  S' ,  l] 

»  -  2  bx  [tanh  ( +5’)]  (2.15) 


which  is  the  result  given  by  Langmuir  for  the  semi-infinite  plasma  bounded 
by  an  emitting  plane. 

Thus  far,  we  have  been  concerned  with  electron  rich  cases.  The  ion 

rich  cases  can  be  derived  simply  be  re-defining  the  parameters  a,  g, 

ti  and  k  as  shown  in  the  table  below. 

LKJ 


Symbol  Electron  Rich 

g2  n+(L)/n_(L) 

a2  n+(0)/n_(0) 

XD0  KT  V^0)) 

T|  eV/KT 


Ion  Rich 
n_(L)/n+(L) 
n_(0)/n+(0) 

*T  «c/(e%(0)) 
-  eVAT 
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With  the  new  definitions,  the  Poisson's  equation  becomes  identical  to 
Eq.  (2.2)  with  a  <  1. 

B.  NUMERICAL  EXAMPLES  WITH  PURE  METAL  EMITTERS 

To  give  some  idea  of  the  types  of  densities  and  potential  profiles 
which  can  be  expected  inside  a  "black  body,"  we  give  some  numerical  re¬ 
sults  based  on  emission  data  of  pure  refractory  metals.  In  practice 
the  quantity  p  =  n+(L)/n  (L)  is  not  easily  obtained  because  the  data 
on  ion  emission  rate  from  various  metals  are  scarce.  For  the  numerical 

computations  we  use  the  following  empirical  formula  on  ion  emission  from 

29 

molybdenum  emitter  as  given  by  Wright  '  and  the  Richardson  equation  for 

2 

electron  emission.  Current  densities  are  given  in  amps  per  cm  . 

e0+  k 

in  J+  =  28.39  -  —  -  0.453  in  T  +  6.22  X  10  T  ,  (2.l6a) 

J_  =  AT2  exp(~  (2.l6b) 


with 

2 

A  =  55  amp/cm 
0  =  8.3  volts 

0_  =  4.37  volts  . 


In  the  computation  the  emitters  are  assumed  to  be  1  cm  apart.  In  Fig. 
2.4  plots  of  n  ,  n  (0),  n  (0)  and  X  versus  1/T  are  given.  The 

O  t  “  DU 


30 


The 


neutral  density  nQ  is  obtained  from  the  data  provided  by  Honig.- 
graphs  rhow  that  above  2000  °K  the  space  charge  is  practically  neutral¬ 
ized  at  x  a  0.  For  T  <  1800  °K  the  diode  space  consists  almost  en¬ 
tirely  oi  electrons,  the  density  is  very  low  and  the  Debye  length  is  of 
the  order  of  the  diode  separation.  Figs.  2.5  and  2.6  give  potential  and 
density  profiles  for  T  =  2200  °K. 

We  were  somewhat  startled  to  find  that  the  space  charge  neutralization 
requires  only  a  small  amount  of  ion  emission  from  the  walls.  For  the 
example  shown  in  Fig.  2.6,  the  ion  density  at  the  wall  is  10^/cc,  which 


FIG.  2.4.  ELECTRON  DENSITY,  n_(0),  ION  DENSITY,  1^.(0), 
DEBYE  LENGTH  Xqq  IN  MID-PLANE,  AND  THE  NEUTRAL  DENSITY 
n„  AS  A  FUNCTION  OF  1/T  FOR  A  PURE  MOLYBDENUM  EMITTER 
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FIG.  2.%  ELECTROSTATIC  PCTXJfTLAL  ts  =  eV/*T  AS 
A  FC3CTI0*  OF  THE  DISTANCE  FHOW  W^-PULXE. 


-G  2 

corresponds  to  an  ion  current  density  of  10  icps/cn  .  Even  with  such 
a  snail  ion  emission  rcte,  we  can  derive  a  neutral  plana  whose  density 
is  orders  of  nagnitude  greater  than  that  predicted  by  the  analysis  of 
Nottingham  (ions  absent)  or  Eicbenbats  and  Hemqviet  (nonequiiibriua 
analysis  with  unfilled  regions  in  phase  space) . 

The  midplane  electron  density,  ion  density  and  neutral  density  versus 
i/T  for  tungsten  and  niobium  are  shown  in  Figs.  2.7  and  2.8.  For  tung¬ 
sten  the  data  on  the  neutral  density  are  obtained  froa  graphs  provided 
by  Honig,  end  the  following  eupirical  formula  as  given  by  Saith^1  is 
used  to  caapute  the  ion  density  at  the  tungsten  emitter. 


o 


aw  nc*cc~'  »v  izrri 

■jt  t.R  "tit  (Si  "Ttlu  ICW'JT; 

£T  T13  .«  1CTTO  "TT 

fotaca  wcTzma  fa^^rra  .  c* 


-leal 


For  tfaa  |eoenl  case  of  nonzero  3,  on*  nut  resort  to 
integration.  Ibe  following  nonaallsrtlon  Is  found  to  bo  convenient  for 
nunerical  analysis 


y*T  c 

TT 

e  a  is; 


The  variable  P"  asitnrts  the  radius  ia  units  of  ^//j,  uhleh  is  a 
kacara  quantity,  ahile  ths  variable  P  is  'normalized  to  1^,  nblcfc  is 
aoi  kaotc  a  priori.  Poisson's  equation  son  becomes 


d^r  1  dr 

a*-2  *  **  «** 


_1 

a 


jr> 


.e 


-  a^e 


(2.22) 


or  by  re-arraaglag  the  ri|bi  hasd  side  *e  have 


-L  _i_  *H\ 

p*  dp*  V  tp'i 


slab  -  ♦  £  e 


(2.22) 


e» 

■O 


«  1 


Bquatioa  (J.2;)  introduce*  a  parameter  ifeiefe  is  generally  a  smell 

zaaSber  directly  proportional  to  the  fractional  donation  fron  neutrality 
at  p*  -  C.  For  oanputa . ioaal  purposes,  tc  used  Bq.  (2.23).  For  » 
gives  nine  of  3,  s  series  solution  esa  be  obtained  for  small  values 
of  .  Let 


V  *  a 


a^f 

£ 


■*  a. 


Tbes 


-‘•a.  ♦  9a.;'  ♦  lea  ; 
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slab  T  <  a.* 


2x*  1 


foe  . 


2  « 


1  •«  •  X 

*»  - 15  *  1  •; 


*t  *  is 


U  -  —£ 


m. 


~S  jia  -  A< 


-  ♦  :fi 


CfU(  Oit  totoUoa  as  tte  simuc  peixt,  tte  fiSfcmmol  aoiatsa  u 
isUfrat«4  anmcalir  toy  *  SfMJCOOL  yncoAn  catl*eaC  '  Ifim  Praetite 
Cormlioe  tnilaklc  at  Ua  S:asjor<  CaquUtiot  Ceiur. 

procedure  Is  repn:«<  at t±  cu. i.eret;  ra^ata  «:  ;  u::l  u« 

bouadorr  ooudiUss  at  tie  sail  is  wu*::o£.  Tie  raeu.lt  tig  x’.rrtji: 
profiles  for  a  aoliatemt  eat  tier  are  sisa-x  stti  tie  topcnisre  as  a 
paraatter  ia  rij.  2.1..  Tie  tftaneter  of  tie  cat t  ter-s  is  tAcs  tc  t« 

1  aa.  lie  car-res  Afaoa  liai  tie  tetct;  astf  potetiei  profiles  Jt-r  ue 
cjUitfrtcal  sysla*  look  like  tic.se  of  tie  plaxar  system  tart  tie  racst 
is  mod  larger  tiax  tie  >*?t  ies^tt.  ' 

2** 
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Incut*  i  ft  the  will  »r'Ke  tove red  by  a  no  no- layer  of  adsorbent. 
Shader  steady  state  conditions,  ion  and  a  too  emission  rates  are  equated 


To  obtats  Ue  steedv  state  toTenf*,  ve  asrune  y  =  and  the  atom 

caitsiM  rate  is  equated  to  the  atom  arrival  rate  which  is  conputed  froa 
be  vapor  pressure  of  the  adsorbent  naterial.  In  this  nanner,  the  steady- 
state  fcce'n*  9  is  determined  and  hence  esission  rates  can  be  conputed 
Strictly  spceving.  <mt  computations  violate  the  thernodrnasi c  equilibrium 
ceaditio*  since  ae  m«B«  that  the  vapor  temperature  is  not  necessarily 
ep sal  to  the  rrtal  reapers  tart.  Ctuy  the  resultiz^  emission  data,  ve 
ccapste  the  4r_<ty  and  potential  profiles  for  several  differed  sets  of 
ccaatiaa.  Stoat  u  Fiji.  2.11  and  2.L2  are  the  profiles  for  a  typical 
c«m  at tr  cesian  vapor  and  leaptes  emitters.  Sole  that  the  space  charge 
is  neutral  i*  the  greeter  portion  of  the  diode  space  despite  the  large 
difference  sx  the  electron  and  aox  densities  at  the  vail. 


fsC-.  2.11.  TESCKETICtL  SC  POTEyiUl.  PROFILE  FOP.  A  DIOEE 
1PC5IP  VALLS  A53>  cssrcx  VAPC6 .  For  the  coneputa- 
t  ica.  the  vail  tempera  tare  vas  1  ~<S.  “1  and  the  vapor 
tenperaiare  vas  bl  °C.  The  vail  separation  vas  takes 

to  he  1  a*. 


FIG.  2.12.  ION  AND  ELECTRON  D’SNSITY 
PROFILE.  The  conditions  are  iden¬ 
tical  to  those  of  Fig.  2.11. 


The  midplane  densities  for  sufficiently  large  plate  separation 
(L  »  \  )  are  computed  for  tantalum  emitters  with  sodium,  potassium 

and  cesium  vapors.  Shown  in  Figs.  2.13,  2.14  and  2.15  are  the  midplane 
plasma  frequency  and  the  Debye  lengths  as  functions  of  emitter  tempera¬ 
ture  with  the  vapor  temperature  as  a  parameter.  In  the  range  of  emitter 
temperatures  covered,  the  Ta-Na  system  is  electron  rich,  the  Ta-K  system 
is  either  electron  or  ion  rich  depending  on  the  emitter  temperature,  and 
the  Ta-Cs  system  is  ion  rich.  It  is  interesting  to  note  that  even  with 
sodium,  whose  ionization  potential  is  higher  than  the  work  function  of 
tantalum,  a  good  workable  plasma  can  be  obtained.  The  probability  of  a 
sodium  atom  being  ionized  by  tantalum  is  small,  but  a  sufficiently  high 
number  become  ionized  to  neutralize  'he  electron  space  charge  in  the 
diode  space. 
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2.13.  PLASMA  FBBQDEMCY  AXD  THE  QB8TK  LE3K7IH  FOB  A  SOCfTOM 
LASMA  PRODUCED  BETTSES  TTO  TAXTALGM  PLATES  VS  PLATE  IS0DU- 
rcaE.  S(MtliB  temper* tore  Is  used  as  a  parameter. 


FIG.  2,1'-.  PLASMA  FKEQCE9CT  AXD  TEE  DSBTE  LZ3GTB  FOB  A 
POTASSIUM  PLASMA  PRODOCED  BETTIES  TTO  TAXI  ALEX  PLATSS . 
The  quantities  are  plotted  as  Tuoctloas  ol  taataloa 
temperature  trlth  potassium  vapor  temperature  as  a 
parameter. 
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2W'fs?*  TO*E«UC  ** 

TIG.  2.K.  PLASM  HmCI  JUSD  TO  OMETt  LXMn  RK 

a  cssrar  Tusau  womo  »p»m  no  wtuct  runs. 

The  qsiaUUM  are  fires  as  Sxscllsas  of  uauiai 
taoentare  el  :i  teslas  vapor  twytJilert  as  a 
parameter. 


C.  005 m* 


vm  senna  «usu 


Sodiaa  pi 


is  prodaced  is  a  dtlode  sfclct  is  Ocscrlted  la  tie 


sdeuilc  fracing*  of  fip.  2. It  sad  2.1*.  Tie  diode  consists  of  iw> 
circular  ustalas  sot  loss  of  diameter  1 laches  separated  by  1  tcch. 
The  batten  are  heated  by  electros  bujsrlKat,  sad  the  stole  tube  is 
Unersed  it  as  oil  bath  to  keep  the  mb?  tempera  tore  uni  lent  sad  loa. 


Prorisios  is  eede  sin  to  vary  the  oil  betb  t 


tare.  The  teo  giidded 

.1  -  -----  3-  TV> 


probes  are  for  the  Uafaa  ssre  meesur  aments  discussed  el  sect:  ere.  The 

density  is  measured  by  s  Uspdr  probe  consisting  cl  a  long  wire  of 
15/1 ZOZ  inch  diameter.  The  probe  is  provided  with  as  accurate  meter 
shies  indicates  the  probe  distance  freer  the  cesster  line,  sad  provisions 
are  made  to  pereit  angular  notion  of  the  probes  ss  sell.  The  but  toss  are 


surrounded  by  heat  shields  shies  also  act  as  guard  rings. 


.0  confine  the 


plasma  ia  the  radial  direction  and  make  the  diode  closely  approximate  as 
infinite  system,  es  axial  magnetic  field  is  applied. 


J 

1 

\ 
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fig.  2.i£.  scininc  munm  or  m  plasma  ttbe.  a  =  sodtua 

rmpor  Kwrct;  S  s  uatalm  putts;  C  a  {ridded  probe;  D  a 
beet  shield  ({ward  rlap);  E  a  bowbwrder  fllameot;  Vg  » 
boaeerder  power  supply;  V  a  filseent  power  supply. 


FIG.  2.1".  3CHDMTIC  JEASISC 
OF  THE  PHOTO  SYSTEM.  A  a 
tantalus  plate;  B  a  {ridded 
probe;  C  a  Laarcuir  probe. 
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(3.1=) 

*w«  Uroilfcoot  this  section  ve  take  v  >  C.  To  obtain  the  boundary 
conditio*,  ve  shall  as sue*  that  a  certain  fraction  S(v)  of  the  particles 
lapiaciac  «*!»«  the  vail  vith  velocity  v  is  specularly  reflected,  and 
the  rest  perfectly  absorbed  by  the  vail.  In  other  vords,  at  x  =  -  L/2 


fl[P'-  |»v)  =  *(v)f,jp,-  |,-v)  (3.16) 

and  at  z  =  L/2 

fi[p4-*v) =  r(v)£i[p4v)  (3-it) 


39 


uwfflQ 


e*a  f**5  cm  kt  a  capla  »ctr  wfcicfr.  •ilm  um 
M«7  Is  itflfctiss,  9y  ttew  konAsty  to  [}#1^ 

•»  (3-1*)#  cm  mUtiM  f  *p-  *  S^2f  5  «•)  to  efctsia 


*«-*•*•' 


v  ¥]  ■£ 

I  ’  *  -J 


V*  r"3V 


fc  ~  4 


f;  .- ■ 


-feKr.fr-f  -  ■‘•4 

*  '  -{■  •  «*  —j-  ¥]} 


3-19/ 


Fro*  XmeU'i  equation  we  have 


Ij[p)  *  JjCp,*)  +  p«0®[p,*) 


=  P*c®Cp»x)  ~  e  f°  vjf1CP,*,v)  -  *lCp^)}-r  (3-22) 

where  the  total  current  density  I ^  in  the  one  diaensioasi  system  is 
known  to  be  independent  of  x.  The  quantity  I ^  represents  the  driving 
current  which  was  assumed  to  be  zero  In  the  study  of  the  lossless  plasma 
resonance  in  Section  IIIA.  By  using  Eqs.  (3.15)  *nd  (3.19),  £q.  (3-2C) 
can  be  expressed  in  the  following  fora 

/•L/2  ,L/2 

1  »  p«  E[p,x)  +  I  a(lx-5l)E[p,§)d5  +  I  P(x,S)£p,S)<»S  (3-21) 
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tiiid  is  ike  result  obtaiacd  H  iili/ 

Shoes  ;a  Figs.  and  * are  the  reel  and  uupatn  parts  #!  the 
iapediace  for  different  raises  of  5.  Sole  that,  tics  R  =  i,  resoffloace* 
occur  afceoerer  71*1  =  [ZiHlk  2  ncrc  s  is  as  integer .  fawicr,  at  :i 
B  z  C  these  resonances  disappear  except  for  the  loses?  order  «u  .  Stare 
the  spherical  shell  aodel  does  aot  exhibit  La  sties  daapisg.  the  loss  tx 
this  case  can  be  attributed  entirely  to  the  esd  plate  absorption,  xatf 
the  "Q~  of  toe  resonance  is  dtirecilj  related  to  the  reflection  coeffi¬ 
cient  R. 
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FIG-  .  HBKALIZZX)  UffCnuaCE  PD®  A  PAIR  OF  IXFIXITE 

Gains  npnasED  is  as  ixfiwte  plasma. 


A  *« 


:ke  Co*  mm*sc«  ix  rwncu.  First,  t*  gri  £  system  exhibits 
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c.  comojxist  awe 

Electrostatic  nsoaaaca  ix  a  aoscuioro  sitsu  Codt  are  predicted 
in  its  lag  tit  hydrodyxaouc  model  aad  a  ficlUioes  ixwnfttn  coaCnoe. 

Tbo*  resoaaacas  pobtii/  will  be  modified  greatly  ties  Atactic  models 
are  used  to  study  tie  diode.  lit  asalysis  of  a  aaifcrt  piata  diode 
otiat  ti«  spherical  shell  velocity  distribution  shoes  that  wall  absorp- 
tloa  plays  aa  important  role  ix  the  de  terminal: on  of  the  diode  impedance. 
Ic  fact,  shea  complete  ahaorptioa  occurs  at  Use  walls,  the  only  resonance 
of  significance  occurs  when  the  diode  separation  is  exsctly  equal  to  a 
half-wavelength.  Because  of  the  absorption  effect,  resonances  are  not 
found  when  Multiples  of  ha If -wavelength  "fit"  into  the  diode  space. 

Thus,  soee  of  the  resonances  predicted  by  the  hydro- dynamic  wodel  nay 
not  be  observable,  especially  those  found  in  the  electron-rich  diode 
where  all  of  the  electron  trajectories  hit  the  wall.  (See  the  phase 
diagris  of  Fig.  2.2).  However,  the  resonances  predicted  for  the  ion- 
rich  c?se  say  be  observable  since  most  of  the  electrons  in  this  esse 
are  trappod  in  the  diode  jpace  and  do  not  strike  the  wall. 

Neither  the  hydrodynamic  model  nor  the  spherical  shell  model  is 
capable  of  exhibiting  the  Landau  damping  effect.  An  analysis  employing 
*  more  realistic  gaussian  velocity  distribution  will  further  modify  th.' 
resonances  and  the  impedances. 
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The  Boitxaeszs  e^atns  for  the  ;  species  of  the  plane  i*  electric 
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(i/36  =  -)  aatf  uiii form  is  the  direction  of  the  externally  applied  dc 
magnetic  field  -3.  sx  =  :}.  Any  faactiw  of  the  cosstazts  of  tie  notion, 
such  as  the  total  energy,  satisfies  E<.  Ox*  such  solution.  is 


given  by  ' 
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Tr  i 


where  _  si  i.  is  *  rector  aloof  the  aagnetic  field,  r  is  the 
J  *  J  _  _ 

coordinate  rector  perpendicular  to  the  nagcetic  field,  A(r)  is  the 

rector  potential,  aad  V(r)  is  the  scalar  potential.  As  can  be  seec, 

£4.  (--3)  is  given  in  terns  of  the  total  energy  and  the  total  angular 

wmmsvud  which  are  constants  of  notion.  The  electric  and  magnetic  fields 

can  be  derived  f roc  V(r'  and  A{r)  by 

S  =  7  X  A  (7)  ; 

Z  =  V(7)  . 

We  snail  prescribe  the  solution  to  have  Maxwellian  velocity  distribution 
along  the  axis,  and,  hence. 


+  V  -  C  *r)+i(J.Xr)  .  A(r)  -  vfr) 
2  *3  3  “j 


(4.4) 


By  assuming  further  that  the  diamagnetic  field  produced  by  the  plfsma 
rotation  has  negligible  effect  on  the  dc  magnetic  field,  the  vector 
potential  becooes 


A  =  B  X  r 
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and  Eq .  can  be  reauced  to 
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where  __  =  jejBoj/By.  is  the  cyclotron  frequency  of  the  j  species, 
Equation  (~.p)  can  be  rewritten  as 
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Frota  Eq.  (4.6)  we  can  see  that  the  whole  species  rotates  as  a  solid  body 
with  rotational  frequeccy 

When  Eq.  (4.6)  is  specialized  for  electrons  and  ions,  we  have 
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r?  +  e^rAJ.  (t  3) 


The  density  profiles  and  the  current  densities  can  be  obtained  by  per¬ 
forming  the  appropriate  integration  of  Eqs.  (4.7)  and  (4.8). 
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(4.10) 


Je  =  -  e(%  *  r)n,(r)  , 


J .  =  e(o>.  X  r)n  (r) 
i  v  i  i ' 


(4.11) 

(4.12) 


For  the  special  case  of  complete  neutrality,  the  simple  solution 
yields  a  gaussian  density  profile  with  mean  squared  radius 
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(4.13) 


and  the  electrons  and  ions  each  rotate  as  a  solid  body  about  the  axis 

with  frequencies  cx>e  and  respectively.  The  following  restrictions 

must  be  placed  on  oo  and  a),  to  avoid  densities  which  rise  with 

e  i 

increasing  r. 
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C.  SCFfWOSlTIOX  or  DISPUCSD  ELEMESTASi  SOLLTIOXS 

Cbder  tie  assumption  of  complete  neutrality,  a  store  interesting  type 
of  rolctlofa  is  oi> t * > red .  Such  a  solution  is 


V(r-  1)1" 


0).  +  a j.x  . 

i  i  ci 


(4.26) 

wbere  €  is  a  vector  in  the  r,  9  plane.  Equation  (4.2 6)  is  just  the 
convolution  cf  the  eleaentary  solution  and  the  "displacement"  function 
F (§).  Equation  (*..26)  satisfies  the  Boltzmann  equation,  but  there  is  no 
way  to  handle  the  non-neutral  cases  in  representation.  The  self- 

nagnetic  field  will  also  spoil  the  solution.  To  insure  neutrality  the 
distribution  function  for  the  electrons  must  be  of  the  form 
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where  the  rotational  frequency  co is  related  to  u>^  by  Eq.  (4.13^ 

An  example  of  F(§),  which  may  be  applicable  to  a  practical  situation, 


is 


F(t)  =1  for  0  <  |1|  <  bc 

=  0  for  H|  >  b 


(4.28) 


This  is  expected  to  give  a  solution  which  closely  approximates  a  plasma 
produced  between  two  circular  diode  plates  of  radius  bQ  and  confined  by 
an  axial  magnetic  field.  For  ions  the  expressions  for  the  density  profile 
and  the  current  density  become 
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where  JQ,(r)  is  the  9  component  of  the  ion  current  density,  I  and 
y  1  o 

are  zero-  and  first-order  modified  Bessel  functions  of  the  first  kind, 

N  is  the  number  of  particles  per  unit  length,  and  r  is  as  defined  in 
o  m 

Eq.  (4.13).  From  Eqs.  (4.27)  and  (4.28),  the  electron  density  and  cur¬ 
rent  profiles  can  be  obtained  and  are  found  to  be  related  to  N^(r)  and 

J  ( i )  in  the  following  manner 
u  i 


Ne(0  =  ^(r) 


co 

Jee(r)  Jei<r> 


K31) 

(4.32) 
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Equations  (4.29)  a:ld  (4*3°)  can  be  integrated  readily  numerically,  and 
the  results  for  ox  =  -a)  72  are  shown  in  Figs.  4,1,  4.2,  and  4.3.  The 

1  Cl 

rotational  frequency  o>^  is  arbitrary,  but  the  choice  of  o>^  =  -<dc^/2 
gives  a  solution  with  best  confinement  (smallest  mean  radius).  In  the 
figures  the  parameter  varied  for  the  plots  is 
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n. 


8kT, 


m.cn  . 
i  ci 


(4.33) 


The  graphs  show  what  we  expect  of  a  magnetically  confined  plasma  as  the 
magnetic  field  is  varied.  The  density  profile  becomes  flat-topped  as 
the  magnetic  field  is  increased  and,  at  infinite  magnetic  field,  the  pro¬ 
file  becomes  perfectly  flat-topped  with  a  sharp  boundary.  In  general, 
the  density  drop  at  the  plasma  edge  occurs  roughly  within  a  distance  of 

2r  .  The  current  density  is  highest  where  the  density  gradient  is  highest, 
m 

The  ion  drift  velocity  is  of  the  order  of  the  ion  thermal  velocity  at 

the  plasma  edge.  Solving  Eqs.  (4.21)  and  (4.22)  for  o>e  under  the 

condition  of  T  =  T  and  cn.  =  -<jd  ./2,  we  find  to  « o>  /4  =  -co  /2. 

i  e  l  ci  e  ci  l 


FIG.  4.1.  DENSITY  PROFILE.  Normalized  density  is 

plotted  versus  normalized  radius  for  different 

values  of  r  »  JQkT .  / (m  <xr  ) . 

m  v  i  '  i  ci' 
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Tfcos,  the  electros  drift  velocitf  is  also  of  the  order  cf  the  ion  drift 
Telocity,  and,  hence,  negligibly  small  vhes  coapa red  to  the  electron 
thermal  velocity. 

*-  3-  Cutler^  has  derived  an  expression  very  similar  to  Eq.  (-.20) 
bat  alth  as  entirely  different  method.  After  approximation  Cutler 
arbitrarily  sets 


£ 

It  _ 

1  Cl 


This  value  dees  not  satisfy  the  inequality  (-.1-)  for  real  _  and, 
beace,  is  disallomed  in  our  theory. 

D.  EXPEXIMEVTAL  ®ASrt£5ffXIS  Or  DC  SEXSITY  AXD  DRIFT  PROFILES 

Is  the  sodhua  plasma  produced  is  the  manner  described  in  Chapter  HE, 
ae  attempted  to  measure  the  azimuthal  electron  drift  and  the  density  pro¬ 
file.  For  the  measurement  a  modified  Langmuir  probe  is  used.  The  tip 
of  this  pr»  >  is  Shota  in  Fig.  It  cat.  be  rotated  so  that  the  dc 

irmsity  profile  can  be  measured  by  facing  the  probe  surface  toward  the 
buttons,  »n le  plasma  rotation  uay  be  measured  by  facing  the  probe  surface 
is  the  azimuthal  direction.  The  detailed  circuitry  and  the  probe  tech¬ 
nique  used  to  overcame  various  problems  is  discussed  in  Appendix  A. 
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FIG.  LA>Cin;iR  PROBE 


With  the  collecting  surface  facing  the  azimuthal  direction,  the  probe 
curve  is  obtained  first  with  the  magnetic  field  directed  one  way  £  id  then 
with  the  magnetic  field  reversed.  Presumably  the  azimuthal  drift  can  be 
detected  by  noting  the  differences  in  the  probe  current  for  the  two  cases. 
Shown  in  Fig.  4.5  is  a  typical  set  of  probe  curves  in  the  magnetic  field. 
As  the  probe  potential  is  increased  from  some  negative  value,  the  current 
rises  exponentially  and,  at  the  plasma  potential,  breaks  away  from  the 
exponential  rise  to  form  a  "knee."  We  have  plotted  in  Fig.  4.6  the  cur¬ 
rent  at  the  "knee"  versus  radius.  There  is  a  significant  difference  in 
the  probe  current  when  the  magnetic  field  is  reversed,  indicating  that 
an  azimuthal  drift  is  present.  However,  due  to  the  lack  of  an  appropri¬ 
ate  probe  theory,  the  measured  result  cannot  be  quantitatively  related 
to  the  azimuthal  drilt.  A  similar  measurement  for  ion  drift  is  desirable, 
but  considerable  improvement  on  the  probe  circuits  is  required  first 
since  the  ion  currents  are  of  the  order  of  the  leakage  current  or  less. 

At  180  gauss,  the  radius  rm  as  defined  by  Eq,  (4.32)  is  3  cm.  This 

value  is  much  too  high  to  give  current  density  profiles  such  as  those 

of  Fig.  4.6.  Comparison  of  theoretical  density  profile  was  also  made 

with  the  experimental  density  profile  of  a  plasma  produced  in  a  similar 
30 

manner  by  Cutler.  ^  There  again  we  found  that  invariably  the  experimental 


FIG.  4.5.  TYPICAL  SET  OF  PROBE  CURVES.  The  two 
curves  are  obtained  by  keeping  the  probe  in  the 
identical  position  but  by  changing  the  direction 
of  the  magnetic  field, 
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FIG.  4.6.  AZIMUTHAL  PROBE  CURRENT  VS  RADIAL  DISTANCE.  The  two 
curves  are  obtained  with  identical  but  oppositely  directed 
magnetic  fields.  The  parameters  of  the  experimental  plasma 
were:  T  =  1800  °K;  n  =  3  X  10°/ cc;  =  1.7  mm, 

density  profile  had  a  much  sharper  drop  at  the  plasma  edge  than  the  theo¬ 
retical  profile.  The  discrepancy  emphasizes  the  importance  of  the  elec¬ 
tric  field  which  we  have  ignored  in  our  theory.  The  electric  field  may 
be  set  up  in  the  following  manner.  The  primary  electrons  supplied  within 
the  solid  cylinder,  due  to  their  small  Larmor  radii,  have  a  sharper  boundary 
than  the  ions.  Hence,  charge  separation  occurs  and  an  electric  field  is 
set  up  which  pulls  the  electrons  outward  and  ions  inward.  Since  the 
probe  measures  the  electron  current,  it  is  not  surprising  that  the  pro¬ 
file  is  sharper  than  that  predicted  by  a  theory  which  assumes  complete 
neutrality.  Indeed,  measurements  show  that  there  is  an  electric  field 
present,  though  not  concentrated  near  the  plasma  edge  as  might  be  expected. 
Shown  in  Fig.  4.7  is  the  potential  profile  of  the  column.  An  electric 
field  of  approximately  0.25  volts/cm  is  present. 

For  the  study  of  waves  along  a  plasma  column  undertaken  in  Chapter 
VI,  we  shall  assume  that  the  plasma  column  is  uniform  with  a  sharp 
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V.  WAVES  IN  AN  INFINITE  MAGNETOPLASMA 


In  this  chapter  we  investigate  the  waves  in  an  infinite  uniform  hot 
magnetoplasma  keeping  in  mind  our  final  goal  of  obtaining  results  for  a 
finite  cylindrical  plasma  column.  We  a**e  primarily  concerned  with  the 
derivation  of  the  dielectric  tensors  which  describe  the  plasma  behavior. 
We  will  demonstrate  a  new  method  to  separate  the  plasma  conduction  cur¬ 
rent  into  polarization  and  magnetization  currents.  The  separation  is 
found  to  simplify  the  analysis  considerably. 


A.  THE  DIELECTRIC  TENSOR  IN  A  MAXWELLIAN  ELECTRON  PLASMA 

Let  us  consider  the  linearized  collisionless  Boltzmann  equation  for 
the  electrons.  For  the  moment  we  shall  assume  that  the  ions  are  infi¬ 
nitely  heavy.  Then  the  zero  and  first  order  equations  are  given  by 
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(5.2) 

where  =  f^(t,x,v).  Bv  Fourier  transforming  in  x  and  Laplace 
transforming  in  t,  Eq.  (5-2)  is  reduced  to 


(to  -  ik  •  *)*L  -  |  [ V  x  iQ]  --1=1  (E  +  V  X  B)  •  ~  (5.2) 


dv 


dv 


where  now  f^  =  f^[<n,k,v].  We  introduce  the  following  gyro-tensor 


D(q>)  =  n,n  +  (I  -  n,n)cos  cp  -  n  x  I  sin  cp 


(5.^) 


and  its  integral 


L(cp)  =  C  D(tp')<ip 

J0 


=  n,n  cp  +  (I  -  n,n)  sin  cp  +  n  x  l(cos  9  -  l) 


(5.5) 
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where  n  is  a  unit  vector  in  the  direction  of  the  magnetic  field  and 
"i"  is  the  unit  tensor.  The  detailed  derivation  of  the  solution  to 
Eq.  (5*3)  is  given  in  Appendix  B.  The  solution  is  found  to  be 


f>,k,v]  -  -S-  r  expfeg  +  i  L:  "  •  "]  E  •  D  .  ^2  ap  .  (5. 

1  “c  Jo  l“c  “c  J  8v 


6) 


This  solution  can  be  verified  by  substituting  it  directly  into  Eq.  (5.3). 

44  1  < 

The  above  solution  agrees  with  those  given  by  Fried  and  Stix.  J 
Through  the  use  of  the  gyrotensors  D(cp)  and  L(co),  our  solution  appears 
in  a  nuch  simpler  form  than  that  given  by  Stix.  The  current  density  can 
be  obtained  by  the  Integration 


<  pv  >  =  -  e 


J vf^ojjk/vjd^v  . 


(5.7) 


In  Appendix  C  we  show  that  it  is  possible  to  separate  the  current  density 
into  polarization  and  magnetization  currents,  as  suggested  by  Lorentz's 
electron  theory  of  matter,  i.e. 


<  pv  >  =  icn  P  -  ik  x  M 


(5.8) 


where  D  =  e  E  +  P  and  B  =  u  (H  +  M) .  For  a  Maxwellian  plasma  we 
o  o 

obtain  (see  Appendices  D  and  E), 
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2  Jo 
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^  f*  •  rf 


OJ 
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L(-cp)  .  E  dcp  (5.9) 
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M  = 


tov0 
3  <2 

OJ  c  U 
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B  exp 


iojcp 

co 


4  (k 


So 


dcp 


(5.10) 


where 


*ji(y)  =  (I  -  n,n)cp  sin  cp  +  (n  x  l)cp(l  -  cos  cp)  +  2n,n(l  -  cos  cp) 
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c  «l+4  f*  ;  e-U  dc 
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cos  5-e  dtp 


(5.i3) 


where 


u  =  ^  -  -4  (k  •  L); 
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k  and  k  era  t*»«  m*  sabers  perpeadlcolsr  sad  parallel  to  the 
*  l 

ftti  r  field  respectively.  Xc  Appendix  D  we  skew  bow  to  express  Eqs. 
('.11}  to  {5.13}  la  terms  01  the  Hilbert  transform  of  the  gaussian  as 
deflate  by  Fried  sad  Conte.  ''  Use  result  Is 


rbere  I  =  I  (X*V2)  ere  modified  Bessel  functions. 


a>  -  nui 


c„  = 


n  ft 


k  vc 


k  = 


J2  V 


u> 


and  Zn  =  Z(Cn)  is  the  Hilbert  transform  of  the  gaussian,  commonly 
called  the  Fried  function,  defined  by 


2 

1  r  ~x 

2(cn)=-prJ  ~rr  Im(cn)  >  0  .  (5.17) 

Properties  of  this  function  are  discussed  in  References  3  and  4,  and  a 

few  simple  relations  involving  Z(£)  are  given  in  Appendix  D. 

The  dielectric  tensor  described  by  Eqs.  (5.H)  to  (5.16)  is  much 

simpler  than  the  "dielectric"  tensor  derived  by  earlier  authors  and  de- 

13 

scribed  in  the  book  by  Stix,  J  In  the  representation  given  by  Stix,  the 
plasma  conduction  current  is  not  divided  into  magnetization  and  polariza¬ 
tion  currents  and  hence  the  "dielectric"  tensor  is  not  truly  the  tensor 
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which  relates  the  electric  field  to  the  displacement  field.  Even  under 
quasi-static  approximation,  our  results  do  not  agree  with  that  given  by 
Stix. 


B.  THE  DIELECTRIC  TENSOR  IN  PLASMAS  WITH  TEMPERATURE  ANISTROPY 

The  derivation  of  the  dielectric  tensors  for  plasmas  with  temperature 
anistropy  is  given  in  Appendix  C.  Equations  (5.1*0,  (5.15),  and  (5.16) 
remain  unchanged,  but  some  of  the  variables  must  be  redefined  as  follows 
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The  limiting  case  of  interest  is  Vq^  =  0.  In  this  limit  the  dielec¬ 
tric  tensor  elements  become 
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These  agree  with  the  results  obtained  by  Lichtenberg  and  Jayson.  Lich- 
tenberg  and  Jayson  derived  the  tensor  by  assuming  that  the  plasma  con¬ 
sisted  of  a  continuum  of  streaming  electrons  with  a  gaussian  velocity 
distribution  along  the  magnetic  field.  They  obtained  the  polarization 
tensor  by  integrating  over  the  contribution  from  the  continuum  of  electron 
streams.  In  Chapter  VI  the  dispersion  relation  for  the  Lichtenberg- 
Jayson  model  is  compared  with  that  of  a  fully  thermal  plasma  model. 

The  cold  plasma  dielectric  tensor  is  recovered  by  letting  vQ 

”  II 

approach  zero  in  Eqs.  (5.18)  to  (5.20). 

C.  DIELECTRIC  TENSOR  FOR  A  MULTI -COMPONENT  PLASMA 

The  polarization  vector  for  a  multi -component  plasma  is  given  by 
the  sum  of  the  polarization  in  each  of  the  components.  Thus,  the 
dielectric  tensor  is  given  by 


e  -  i  + 


I  ”j 


th 

where  n .  is  the  electric  susceptibility  tensor  for  the  j  species. 

J  , 

The  elements  of  n .  are  given  by 


2 
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n  .  = - ^  c  .  exp 

ZZJ  a)2  °J 


[4]Uk 


(5.21) 


it  .  =  -it 

xyj  yxj 


io)2  C  • 
Pj  QJ 

0)  .01 

CJ 


exp 


[4]H- 


Vi+Z„-1> 


(5.22) 


and 


„  “4C°J 

xxj  yyj 


exp 


X2"* 

J 


I  Wl  -  Vl>  (5-23) 


-  65  - 


where 


I  = 
n 
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z  =  2(C  J  , 

n  v*nj  ’ 


CD  -  UD 
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Cj 


^  "  /Tk 


II 6J 
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D.  THE  QUASI-STATIC  DISPERSION  RELATION  FOR  THE  INFINITE  PLASMA 

Having  obtained  the  dielectric  tensors,  we  can  now  derive  the  dis¬ 
persion  relation  for  the  plane  electrostatic  waves  in  the  magnetoplasma . 
The  distinguishing  feature  of  electrostatic  waves  is  that  the  electric 
field  is  derivable  from  a  scalar  potential 

E^t)  1.-7  0(x,t)  ,  (5.2*0 

and,  hence, 

V  x  E^(x, t)  =  0  . 

This  approximation,  commonly  called  the  quasi-static  approximation,  is 
valid  whenever  the  phase  velocities  of  the  waves  are  much  less  than  the 
speed  of  light.  Using  this  analysis  and  assuming  plane  wave  propagation 
of  the  form  exp[icDt  -  ik  •  x],  we  obtain  the  following  result,  relating 
the  potential  and  the  electric  field, 

E[cD,k]  =  ik  0[cD,k]  .  (5.25) 

From  the  second  Maxwell  equation  we  get 
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where  k“  *  k^  +  k^  is  iki  iqnrt  of  tbe  eat*  saner  jeryndccUr  to 
the  nagoetic  field.  As  expected,  tie  dispersion  relation  is  independent 
of  the  ancle  beraeea  the  perpendicular  caopdeai  of  k  and  ike  x  am. 

This  dispersion  relation,  ~lo^l  It  appears  to  be  different,  is  identical 

1  - 

to  the  result  obtained  by  Stix.  * 

In  an  electron  plasaca  the  dispersion  relation  can  be  pns  :*  terns 
of  one  integral  or  one  senes  as  ioUon 


are  as  defined  in  Section  YA . 
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aa  tkt  <iipmioi  ralatioa.  ly  usiag  Hj.  (5^9)  and  the  asyaptotic 
mixtion  for  Z  (lrn  in  Appendix  D,  00  obtain  the  follovlcg  dispersion 
relation  for  the  *lmstcla'  asides. 
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This  htftmcft  predicts  raeonaecea  seer  the  ha  mo  ales  of  the  cyclotron 
frequency.  Semrml  papers  on  the  expert nantal  obserrations  of  such 
raaonaatar  have  been  pnhllshed.  A  sort*  net  hod  of  sianaing  the  series 

of  (5.3I;  la  texaa  of  Xielsoe’s  cylinder  function  has  been  given  by 
twflsr.  A  mpster  piugran  based  on  this  result  is  used  In  obtaining 
■Jw?  electric  sneceptiUllty  tensor  for  the  ions. 
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In  the  cylindrical  coordinate  system  we  let 


k  =  k  cos  cp 
x  X 

k  =  k  sin  cp 

y  i 

x  =  r  cos  0 
y  =  r  sin  0 

where  cp  is  the  angle  between  the  r  axis  and  the  plane  containing  the 
z  axis  and  k.  The  potential  can  now  be  expressed  in  the  following 
manner 

0[cn,r,0,z)  =  0[cn,k  ,k^,cp]exp[-ik^z-ik  r  cos  (0-cp)  ]  t  (5.35) 

The  cylindrical  wave  is  synthesized  by  taking  these  waves  with  identical 

amplitude,  k  and  k  but  with  phase  variation  of  the  form 
x  II 

eto,  W*]  ■  el“p 

50 

where  m  is  an  integer.  Synthesis  of  such  plane  waves  yields 


0  fQ+n  ± 

0[cn,r,0,z)  =  —  I  e  “*exp[-ik  z-ik  r  cos(G-cp)  ]dcp 
2*  JQ_n  II  1 


-ik  z 

=  0O  ®  m  6  W*  (5.36) 


where  J  (k  r)  is  the  Bessel  function  of  the  order  m.  Similar  inte- 

m  x 

grations  for  the  electric  and  displacement  fields  yield 
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(5. 40) 
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This  technique  of  cylindrical  wave  synthesis  has  been  employed  by  others 

♦4 

in  the  past.  In  our  representation  e  is  independent  of  <p,  the  inte¬ 
gration  variable.  In  Stix’  representation  che  dielectric  tensor  is  a 
function  of  cp,  and  the  integrals  for  the  displacement  vector  is  some¬ 
what  more  difficult  to  evaluate  then  in  our  case. 

The  solution  obtained  above  is  valid  for  an  infinite  plasma.  Our 
ultimate  goal  is  to  find  the  solution  for  a  finite  plasma,  particularly 
a  cylindrical  column  of  plasma.  In  a  finite  plasma  the  most  difficult 
task  is  the  treatment  of  the  boundary.  To  make  the  problem  manageable, 
we  must  assume  that  the  plasma  column  is  uniform  with  a  sharp  boundary. 
In  reality,  a  sharp  boundary  cannot  occur  because  of  thermal  motion, 
though  such  a  model  may  approximate  the  plasma  column,  as  suggested  by 
the  theory  and,  in  particular,  the  measurements  discussed  in  Chapter  IV. 
In  Chapter  VI  we  shall  use  the  fields  given  by  Eqs,  (5.37)  to  (5.42)  as 
the  solution  within  the  cylindrical  plasma  column  and  match  them  to  the 
fields  outside  the  plasma. 
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VI.  ELECTROSTATIC  WAVES  IN  A  CYLINDRICAL  COLUMN  OF  MAGNETOPLASMA 


In  this  chapter  we  investigate  the  dispersion  relation  for  electro¬ 
static  waves  propagating  along  a  column  of  cylindrical,  collisionless 
magnetoplasma.  The  geometry  we  consider  is  that  of  an  infinitely  long, 
uniform,  cylindrical  plasma  enclosed  within  a  cylindrical  waveguide  as 
shown  in  Fig.  1.1.  First,  the  dispersion  relation  for  the  electrostatic 
waves  is  derived.  Some  numerical  solutions  to  the  dispersion  relation 
for  a  Maxwellian  electron  plasma  are  presented  and  compared  with  the 
dispersion  relations  based  on  other  plasma  models.  Our  results  are  also 
compared  with  the  available  experimental  dispersion  relations.  At  low 
magnetic  fields,  the  surface  wave  is  found  to  be  unstable,  and  the  insta¬ 
bility  character  is  investigated.  The  unstable  ion  surface  wave  is  also 
studied. 

A.  QUASI-STATIC  DISPERSION  RELATION 

Since  there  are  two  distinct  regions  within  the  waveguide,  we  need 
two  sets  of  solutions,  one  valid  within  the  plasma  and  the  other  valid 
in  the  free  space  region  between  the  plasma  and  the  conducting  wall.  For 
the  plasma  region  we  shall  use  the  quasi-static  solution,  obtained  in 


the  previous  chapter  and  given  by 

Ez[<u>r>6>z)  =  “ik  0O  exp[im  6  -  ik(z]Jm(k^r)  ,  (5*37) 

im  0 

E0C<u,r,9,z)  =  r  °  exp[im  9  -  ik^jj^k^r)  ,  (5-33) 

E  [o),r,9,2)  =  k  0  exp[im  9  -  ik  z]-J’(k  r)  (5-39) 

r  i  o  II  m  1 

D2C<u>r>e>z)  =  -ik1(e2Z  0O  exp[im  9  -  ikz]Jm(k^r)  ,  (5. to) 

rim  e  *| 


D  [<jj,r,9,z)  =  0  exp[im  9  -  ik  zj  J  (k  r)  +  k  e  J'(k  r)  ,(5.4l) 

r ’  ’  ’  '  o  II  Jl  r  mv  x  '  x  xx  m'  x  ' 

rim  e  -j 

De[<u,r,9,z)  =  0q  exp[im  9  -  ik(z]  — ^  Jj^r)  -  k^e^J^k^r)  (5.42) 
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In  the  free  space  region,  b  <  r  <  s,  the  solution  to  the  equation 

p 

V  0  =  0  with  the  requirement  of  0(a,8,z)  *  0  is  given  by 

lra6-ik  z 

0(r,8,z)  a  A  [j  (ik  r)H  (ik  s)  -  J  (ik  a)H  (ik  r)]e  11  (6.1) 

'  ’  *  7  o  m'  ||  7  m'  ]|  '  m'  ||  7  m'  ||  7  J  '  7 


where  J  (ik  r)  and  H  (ik  r)  are  Bessel  and  Hankel  functions  of  order 
m.  The  components  of  the  electric  field  are  then  obtained  from  Eq.  (6.1) 
by  E  =  -70.  At  the  free  space  plasma  boundary,  the  tangential  component 
of  the  electric  field  and  the  normal  component  of  the  displacement  vector 
are  matched.  This  yields  the  following  dispersion  relation  for  the  waves 


and  the  dielectric  tensor  elements  are  as  described  in  Chapter  VC.  The 

dispersion  relation  of  Eq.  (6.2),  in  appearance,  is  very  familiar.  In 

7  q 

the  cold  plasma  model  and  in  the  Lichtenberg-Jayson  (T  =0)  model,7 

the  dispersion  relation  is  identical  to  Eq.  (6.2).  The  only  difference, 

though  by  no  means  small,  is  in  the  dielectric  tensor  elements. 

To  solve  the  dispersion  equation  in  the  form  given  by  Eq.  (6.2)  is 

a  formidable  task.  For  example,  the  computation  of  the  dielectric  tensors 

above  involves  infinite  summations  of  products  of  the  Bessel  and  Fried 

functions.  Only  under  certain  conditions  do  these  summations  converge 

rapidly  enough  for  computer  analysis.  In  restricting  our  study  to  the 

cylindrically  symmetric  modes  (m=0) ,  we  have 


(6.3) 


Even  for  this  case,  a  thorough  analysis  of  all  the  modes  is  formidable. 
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B.  THE  CASS  Of  HIGH  MAGNETIC  FIELD  (co  >  cq  )  WITH  STATIONARY  IONS 

v  ce  pe' 

1 .  Comparison  of  the  Pj«pcr»< on  Eolations  for  Several  Different 
Plasma  Modal a 

The  dispersion  relation  for  the  electroetetlc  saves  In  a  cold 

7 

pi a ana  column  saa  first  derived  by  Trlvelpiece  and  Its  solutions  are 
cosaonly  called  the  Trlvelpiece  sodes.  The  hydrodynamic  model,  obtained 

by  taking  the  momenta  of  the  Boltzmann  equation,  saa  used  by  Agdur  and 

3 

Veiasglas  to  study  the  temperature  effects  on  the  Trlvelpiece  modes. 
Unfortunately,  the  hydrodynamic  model  does  not  account  properly  for  the 
temperature  effects.  In  this  section  se  present  the  dispersion  relation 
for  saves  in  a  Maxselllan  plasma  column  and  compare  it  sith  the  disper¬ 
sion  relations  derived  by  using  other  models.  Ve  shall  assume  that  the 
ions  are  infinitely  heavy. 

Equation  '6.3)  is  solved  numerically  for  the  Maxselllan  electron 
plasma  by  as aiming  real  co  and  seeking  complex  k^.  The  computation 
procedure  is  described  ir  Appendix  E.  The  result  for  the  lowest  order 
mode  is  shown  in  Fig.  6.1.  Also  shown  on  the  same  graph  are  the  disper¬ 
sion  relations  for  saves  in  a  cold  plasma  colunm,  in  a  ’’hydrodynamic” 
plasma,  and  in  an  infinite  Maxselllan  plasma.  For  the  infinite  hot  plasma, 

the  dispersion  relation  presented  ie  the  lowest  order  Landau  node  a a 

42 

defined  by  Sinonen. 

At  los  values  of  a),  the  three  models  yield  identical  dispersion 

relations,  but,  at  high  values  of  u>,  the  dispersion  relation  for  the 

Maxselllan  plasma  column  merges  smoothly  into  the  hot  infinite  plasma 

dispersion  relation,  shile  the  dispersion  relation  for  the  hydrodynamic 

model  approaches  the  thermal  speed  of  the  one  dixencional  adiabatic  Bohn 
51 

ana  Gross  saves.  As  se  can  see,  the  save  in  the  column  of  the  Maxsel¬ 
llan  plasma  is  Landau  damped. 

The  dispersion  relation  for  the  zero  transverse  temperature 
model  (T  =  0)  was  also  computed,  but,  at  the  high  magnetic  field  con¬ 
sidered  above,  the  results  for  T  ;  0  and  T  =*  T  are  identical. 

i  i  li 

This  is  not  surprising  since,  at  high  magnetic  fields,  tne  transverse 
motion  of  the  electrons  is  Inhibited. 
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Their  experiment  was  conducted  in  a  duoplasmatron  type  hydrogen  arc 
source,  which  produced  a  plasma  with  high  longitudinal  temperature.  Be¬ 
cause  of  the  high  longitudinal  temperature,  the  model  with  =  0  is 
well-suited  for  the  analysis.  The  experimental  density  profile  given 
by  Maiuberg  end  Wharton  is  highly  nonuniform  as  shown  in  Fig.  6.2.  For 
computational  purposes  the  plasma  is  approximated  by  a  uniform,  cylindri¬ 
cal  plasma  column  as  shown  in  curves  A  and  B  in  Fig.  6.2.  For  the  uni- 

Q  o 

form  column,  the  density  is  assumed  to  be  1.6  X  ICT/cm  ,  which  matches 
the  asy^>totic  behavior  of  the  theoretical  result  to  that  of  the  experi¬ 
mental  curve.  The  radius  for  column  A  is  chosen  so  that  the  total  number 
of  electrons  in  the  column  equals  the  total  number  in  the  nonuniform 
plasma.  Since  Malmberg  and  Wharton  report  that  the  absolute  value  of 


FIG.  6.2.  DENSITY  PROFILE  FOR  KALXBERG  AND  WHARTON  EXPERIMENT. 


the  density  profile  may  not  be  accurate,  the  procedure  used  to  obtain 
the  equivalent  radius  may  not  be  valid.  In  fact,  better  agreement  with 
the  experimental  dispersion  relation  is  obtained  when  the  column  is  as¬ 
sumed  to  be  2.25  cm  in  radius  as  shown  in  curve  B.  From  the  data  supplied 

by  Halmberg  and  Wharton,  we  find  =  0.1Q  cm  snd  ca  /(u  =  4.0.  The 

D  C  p 

real  part  of  k^  obtained  via  the  dispersion  relation  is  compared  with 
the  experimental  result  in  Fig.  6.3.  The  curves  marked  A'  and  B'  are 


A'  «  plasma  radius  b  =  3-5  cm 
B'  =  plasma  radius  b  =  2.25  cm 

FIG.  6.3.  ANGULAR  FREQUENCY  VERSUS  WAVE  NUMBER. 
Experimental  data  obtained  by  M8lmberg  and 
Wharton  is  compared  with  theoretical  results. 


the  theoretical  results  for  plasmas  represented  by  curves  A  and  B  of 
Fig.  6.2.  In  Fig.  6.4  we  have  a  graph  of  Im(k  )/Re(k^)  versus  v^/vq, 

where  v  =  a/Re(k  )  is  the  phsse  velocity.  The  theoretical  curve  is 

P  i! 

computed  by  using  B  of  Fig,  6.2.  The  agreement  between  the  theory 
and  experimental  points  is  good,  despite  the  approximate  nature  of  the 
model . 
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FIG.  6.4.  kt/kr  VERSUS  (vp/v0)2  FOR  MAIM3ERG  AND 
WHARTON  EXPERIMENT. 

C.  SURFACE  WAVE  INSTABILITY  IN  ELECTRON  PLASMAS  AT  LOW  MAGNETIC  FIELDS 

When  a)  <  cd  ,  the  lowest  order  mode  in  s  cold  electron  plasma 
c  p 

column  is  a  surface  wave  in  the  range  cdc  <  cd  <  o^,  where  cd^  = 

\J (cd2+  CD2)/2  (see  Fig.  1.2).  When  the  dispersion  relation  [Eq.  (6.2)] 
was  first  solved  in  this  region  for  real  cd,  instability  was  indicated. 
In  this  section  we  establish  the  presence  of  unstable  surface  waves  and 
study  the  character  of  the  instability. 

One  way  to  determine  the  presence  of  an  instability  is  to  map  the 
negative  real  k  axis  into  the  complex  cd  plane  and  see  if  the  roots 
fall  into  the  lower  half  plane.  Such  a  mapping  for  a  typical  surface 
mode  is  shown  in  F eg .  6.5.  Instability  is  indicated  since  the  roots 
fall  below  the  real  cd  axis  for  a  range  of  kb. 
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FIG.  6  3.  NEGATIVE  REAL  k..b  AXIS  HAPPED  INTO  THE  COMPLEX 
cu/cup  PLANE  VIA  THE  DISPERSION  RELATION,  SHOVING  UNSTABLE 
ROOT. 

All  Instabilities  can  be  classified  into  two  types,  convective  and 
nonconvective  (absolute)  instabilities.  A  nonconvective  instability 
grows  with  time  until  limited  by  nonlinear  effects,  while  the  convective 
instability  is  one  in  which  a  perturbation  is  amplified  in  space.  We 
have  established  that  an  instability  is  present,  but,  as  yet,  we  have 

not  determined  whether  it  is  convective  or  nonconvective. 

14  15 

Derfler  ’  developed  two  methods  of  classifying  instabilities.  The 
first  method  is  to  map  the  complex  cn  plane  into  the  complex  k  plane 
through  the  "positive  frequency  part"  of  the  dispersion  relation  and  study 
the  saddle-points  (cbo/dk  =  0),  which  are  caused  by  the  merging  of  the 
two  roots  of  the  dispersion  relation.  The  way  in  which  the  roots  verge, 
as  the  imaginary  part  of  the  frequency  is  changed,  is  intimately  related 
to  the  classification  of  the  instability.  Nonconvective  instability 
occurs  if  and  only  if  the  two  roots  merge  across  the  negative  real  k 

s; 

axis  with  frequencies  in  the  lower  half  plane.  This  criterion  is  more 
difficult  to  use  than  the  second  one  which  is  described  below. 

The  saddle-points  when  mapped  into  the  complex  go  plane  through  the 
dispersion  relation  become  branch  points.  If  the  branch  point,  asso¬ 
ciated  with  the  instability,  falls  in  the  lower  half  plane,  the  instability 
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The  .»UMu  .t  lmtlfttei  shewn  la  a  relatively  slid  case  compared 
u  than  at  lever  aaiattic  field*.  Atlanta  vere  amde  to  apply  the  ata- 
Mli it  cn terlae  at  lever  aaprtic  fields,  bet  ee  were  unable  to  locale 
the  tract  pent  doe  to  coa^wtatiomal  diff lenities.  Thus,  whether  the 
mtchallty  resa.xs  oaectte  or  not  for  all  values  of  the  aagnetic 
field  is  still  ax  epee  tmtiw. 

*e  nan  seed?  the  kdtanor  of  the  instability  as  s  function  of  the 
re  name  parameters.  Start  :t  Fi|.  are  napping*  of  the  negative  real 
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sate  mStm.  h  t.  into  the  complex  _  plane  for  the  loses t  order 
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bl>cw  Stti  _ _  __  as  a  parameter.  The  instability  is  strosqj  near  the 

V 

ntlo'-rai  freqseac y  sad  at  -.he  iTtnd  frequency.  Also,  whenever  the 

rrciotros  tanoucs  are  Oeloa  __  =\/  *  „£  2.  the  instability  is 

TV?  c 

enhanced  at  these  tanoucs  as  seer  ic  Fi{.  ;.l.  At  large  negative 


FIG.  6.9.  MAPPING  OF  THE  NEGATIVE  REAL  (k,.b)  AXIS  INTO  THE 

COMPLEX  0)/tUp  PLANE  VIA  THE  DISPERSION  RELATION  FOR 

/cc  =  0.2. 
c  p 


values  of  k^b,  Landau  damping  sets  In  and  makes  the  waves  stable. 

Shown  in  Figs,  6.10  and  6.11  are  the  real  and  imaginary  parts  of  the 

frequencies  versus  k^b  with  0^/00^  as  8  parameter.  Unlike  the  surface 

wave  in  a  cold  plasma  column  for  which  the  dispersion  relation  goes 

asymptotically  to  a)  =  u)y,  the  real  part  of  co  goes  above  as  k^b 

is  increased  indefinitely.  The  magnitudes  of  the  electx'ostatic  potential 

profiles  are  shown  in  Fig,  6.12  for  cjd  /cd  =  0,3.  Note  that  the  fields 

c  p 

concentrate  near  the  edge  of  the  plasma  as  the  frequency  is  raised  above 


The  curves  of  Fig.  6.8  seem  to  indicate  that  the  plasma  is  more  un¬ 
stable  at  lower  magnetic  fields.  However,  the  solution  at  cjjc  =  0  is  known 
to  reduce  to  the  cold  plasma  case  and,  hence,  is  stable.  Figure  6.13 

shows  the  imaginary  part  of  the  most  unstable  frequency  versus  cd  /oo  . 

c  p 

There  are  two  curves,  one  for  the  instability  at  the  cyclotron  frequency 


and  the  other  at  the  hybrid  frequency  o^. 
shown  by  the  dotted  line  was  not  computed. 


The  portion  of  the  curve 
Although  the  instability 


FIG.  6.10.  DISPERSION  DIAGRAMS  FOR  THE  SURFACE  WAVES  SHOWING  THE  REAL 

PART  OF  cd/cd  VERSUS  WAVE  NUMBER  k  ,b. 

P  II 


FIG.  6.11.  DISPERSION  DIAGRAMS  FOR  THE  SURFACE  WAVES  SHOWING  THE 

IMAGINARY  PART  OF  u>/u>  VERSUS  WAVE  NUMBER  k  b. 
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appears  to  increase  as  coc  is  decreased,  the  solution  st  u>c  =  0  is 

known  to  be  stable,  and  thus  a  transition  curve,  such  ss  the  one  shown 

by  the  dotted  line,  must  exist.  With  cyclotron  frequencies  below  a>c  = 

0.1  0)  ,  the  summations  in  the  dielectric  tensor  converge  too  slowly  and 

accuracy  is  lost.  Consequently,  we  were  unable  to  cover  this  range. 

The  effects  of  other  parameters  on  the  instability  are  shown  in 

Figs.  6. 14  and  6.15.  In  Fig.  6.14  we  show  the  effect  of  temperature 

anistropy.  The  roots  of  the  dispersion  relation  are  plotted  with  T  /T 

1  8 

4-003 


s/b>2.22 


ace 


1 

fool 


-* — 003 


FIG.  6.14.  BEAL  WAVE  JOBBER  AS  MAPPED  I  WTO  THE  CCMS>LEX  p  PLAJtE 
VIA  THE  DISPERSION  RELAFICff.  The  locus  diagram  illustrates  the 
•fleet  of  the  transverse  temperature  on  the  surface  wave.  = 

0.5,  a/b  =  2.22  and  =  0.04.  C  ? 


0  01 


FIG.  €.15.  JMPPIKS  or  ns  WBGATIVX  k  b  AXIS  I  WTO  TEE  COMDEX 
a/ffl  pun  WITH  Vg/b  AS  A  PAKAWHS. 


as  •  pmntcr.  Tha  pH— 1  Is  stable  when  T  a  5,  and  the  Instability 

in  its  nonlinear  llalt,  —7  provide  a  — chanl—  to  convert  transverse 

Into  longitudinal  tanperature .  Shown  in  Fig.  6.15  are  tba  roots  of 

dispersion  relation  with  L/b  as  a  para— ter.  Since  we  keep  &-  ant. 

u  ? 
b  constant,  a  change  in  1^/b  is  equivalent  to  a  charge  is  the  overall 

reapers  turn.  We  a—  that  .ut  higher  the  ratio  1^/b  (hence,  Signer 

the  temperature) ,  the  greeter  tbe  instability.  I'tese  graphs  indicate 

tost  the  instability  is  due  to  finite  Lanaor  radii  effects. 

Wo  conclusive  experimental  evidence  of  the  surf  ice  vave  instability 
has  been  found.  There  is,  however,  s  report  of  surface  save  seif  c.tcil- 
lationa  with  aero  magnetic  field  In  a  plsszs  column,' '  Our  tbecr.- 
predicta  strong  instability  si  1cm  Magnetic  field,  but  the  lnstsbili* y 
sbould  disappear  at  zero  magnetic  field.  Since  we  ere  unable  to  compote 
tbe  dispersion  relation  in  the  very  low  ntgoetic  field  region,  «e  csnjsot 
relate  tbe  observed  oscillation  io  our  instability.  Also  at  large  Lanto 
radii,  our  dc  no del  with  sharp  boundary  breads  dcwD. 
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VII.  OOJCLDSICW 


la  our  study  of  electrostatic  wtes  we  caphtsixcd  the  effects  doe  to 
the  finite  boundary  as  well  as  the  effects  of  a*{Wtic  fields  and  plasm* 
temperature .  The  first  part  of  this  report  treated  plane  diodes  is  the 
absence  of  a  magnetic  field.  A  dc  study  of  the  equilibrium  plasma  diode 
vas  undertaken  to  provide  a  basis  for  the  rf  analysis.  The  dc  analysis 
yielded  some  interesting  results  regarding  the  possibility  of  plasma 
formation  within  a  diode.  It  mss  found  tba v  under  equilibrium  coadltico 
plssma  is  formed  eves  when  the  generation  rate  of  one  of  tke  charged 
species  st  the  walls  is  orders  of  magnitude  loser  than  that  of  the  other. 
The  formation  of  plasma  doe  to  the  trapping  of  rare  species  in  the  coulomb 
field  of  the  sbuadsct  species  mss  demonstrated  in  the  case  of  metal  emit* 
ters  where  the  emission  rate  of  the  ions  was  very  muck  loser  thaa  that  of 
the  electrons. 

An  rf  analysis  using  the  hy drodynamic  model  for  the  acmxl font 
plasms  revealed  the  existence  of  a  series  of  electrostatic  re sconces  in 
the  diode.  However,  sc  analysis  of  a  uniform  plasm*  diode  using  a  line- 
tic  plasma  model  stowed  that  the  end  plate  less,  tee  to  absorption  of 
the  electrons  at  the  wails,  mtroteces  a  large  loss  mechanism.  The  loss 
may  make  some  of  the  resonances  predicted  iron,  the  hydrodynamic  model 
unobservable.  Since  this  loss  mechanism  is  related  to  the  process  which 
cscses  end  piste  dlf fusion  la  devices  such  as  the  Q-mechine,  the  rf  mea¬ 
surements  of  the  electrostatic  fields  may  yield  information  regarding  the 
cod  plate  diffusion  coefficient. 

The  second  part  of  this  report  dealt  with  the  s'.jii  of  electrostatic 
wave  propagation  along  a  cylindrical  ooluns  of  Kanelliac  plasma.  First, 
a  dc  study  •'f  the  density  profile  and  azimuthal  asrtEi  it  a  cjiiadricsi 
column  was  conducted.  The  theoretical  analysis  of  the  column  under  the 
assumption  of  complete  neutrality  showed  that  si  high  magnetic  fields 
when  the  10:  Larmor  radii  became  much  smaller  than  the  plasma  radius, 
the  density  drop  was  confined  to  the  small  region  si  the  plasm*  edge. 
Experimental  measurements  indicated  that  the  density  crop  at  the  plasma 
edge  is  in  /net  much  sharper  than  that  predicted  from  cur  tbee-ry  .  Fcr 
the  purpose  of  rf  studies,  we  assumed  that  the  plasms  was  art  fora  with 
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APPENDIX  A.  LCm  DENSITY  PROBE  MEASUREMENT 


Id  this  appendix  we  show  some  of  the  difficulties  encountered  in 
probing  a  low  density  plasma  and  the  ways  in  which  these  problems  are 
overcome.  In  our  sodium  plasma  the  typical  density  is  of  the  order  of 

*7 

2  x  1C  /c c,  which  is  lower  than  most  of  the  laboratory  plasmas.  In  the 
probe  measurements  of  a  low  density  plasma,  the  smallness  of  the  probe 
current  is  the  biggest  problem.  The  probe  size  can  be  increased  to 
increase  the  total  current,  but  beyond  a  certain  size  the  probe  disturbs 
the  plasma  too  much.  Also,  the  probe  should  be  as  small  as  possible  to 
make  local  measurements . 

To  obtain  accurate  probe  measurements  it  is  necessary  to  have  the 
"knee'"  of  the  probe  curve  show  clearly.  This  can  be  attained  by  a  guard- 
ring  arrangement  of  the  probe.  The  above  consideration  is  taken  into 
account  in  the  design  of  the  co-axial  probe  shown  in  Fig.  4.4.  Normally 
the  stem  of  the  guard-ring  probe  is  insulated  so  that  the  current  is  drawn 
only  through  the  probe  surface  and  the  co-planar  guard-ring  surface.  How¬ 
ever,  since  a  film  of  sodium  would  quickly  fora  on  an  insulated  surface, 
causing  leakage  current,  we  decided  not  to  coat  the  stem.  Instead  we 
made  the  whole  probe  as  small  as  possible.  The  size  cf  the  probe  is  such 
that  the  total  current  drawn  to  the  whole  probe  at  the  plasma  potential 
is  less  than  i  percent  of  the  total  electron  current  from  the  emitters 
under  normal  operating  conditions . 

Maneuverability  is  another  desirable  feature  of  a  probe.  The  details 

•  /> 

of  the  probe  drivicg  mechaiism  is  described  by  Simenea.  The  probe  moves 
in  sod  cut  as  well  as  rotates  so  that  the  probe  surface  can  be  made  to 
face  either  in  the  azimuthal  or  sxial  direction.  When  an  axial  magnetic 
field  is  applied,  the  density  can  be  measured  by  facing  the  probe  in  the 
axial  direction,  while  azimuthal  drift  can  be  studied  by  facing  the  probe 
in  the  azimuthal  direction. 

Typical  probe  currents  ire  of  the  order  of  a  few  microamperes. 

This  is  a  level  where  leakagr  currents  from  various  sources  begin  to 
interfere  nth  the  probe  curve,  particularly-  leakage  between  the  guard 
rieg  and  he  probe.  Moreover,  a  dc  amplifier  is  required  to  amplify  the 
voitage  across  the  current  pick-up  resistor  to  drive  the  XY  recorder. 
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A  s climatic  drawing  of  the  probe  circuit  Is  shorn  in  Fig.  A.l.  The 
amplifier  is  s  Kelthley  6CA  electroaster  serplifier,  which  fens  the  desired 
characteristic  of  high  input  impedance,  low  noise  and  drift,  however, 
it  is  necessary  to  keep  the  elect roaster  terminals  near  the  ground  as 
shown  in  Fig.  A.l.  Since  s  terminal  of  the  amplifier  needs  to  be  groueoed, 
we  are  forced  to  use  two  sweep  circuits  insteed  of  one  to  keep  the  guard 
ring  volt sge  equal  to  the  probe  voltage.  The  sweep  circuit  is  ihuam  In 
Fig.  A.2.  Two  such  circuits  with  watched  characteristics  but  insulated 
frou  each  other  were  constructed  sod  pieced  in  a  uetsl  box.  When  s 
aaster  switch  is  turned  on,  the  two  circuits  simultaneously  sweep  out  a 
single  identical  saw-tooth  wave  fore.  The  tine  coast act  of  the  sweep  is 
variable  from  1  to  -  seconds  and  the  amplitude  is  variable  frou  1  to  21 
volts. 


FIG.  A.l.  PBOBX  C1BCCIT 
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b.  soLcnos  or  ns  usttiro  sotrauust  kjdatiox 


fs  this  app wntti  w  obtalc  •  solution  of  tho  hoMrlnd  Boltzuaa 
1— tlo»  for  a  aHfon  plma  is  •  ufsetic  field.  The  eolation  has 
toea  derived  by  others  usiqg  several  different  approaches,  but  ve  present 
a  derivation  here  for  the  sake  of  ccopleteness.  The  analysis  can  be 
generalised  to  any  species  of  the  plasua,  but  to  avoid  extra  subscript 
ee  shall  apply  tae  analysis  to  electrons  only.  We  shall  assuae  that 
the  plaaaa  ia  colliaionlesa  and  that  there  is  no  dc  electric  field.  Then 
the  zeroth  and  the  first  order  Boltzmann  equations  are  given  by 

df  (v) 

.'w  -•  o'  ,  * 

i v  x  8  )  •  — - - =  0  (B.l) 

'  o  o* 
ov 


Let  us  nake  the  following  definition 


—  =05  =  n  05 

m  c  c 


where  n  is  r  unit  vector  in  the  direction  of  the  magnetic  field.  We 
Fourier  analyze  the  linearized  equation  in  space  (im  (k)  =  0)  and 
Laplace  analyze  in  time  (Im  (05)  <0  on  L.I.P.)  to  get 


i(o5  -  k  •  v)  f  [oo,k,v]  4  ci5  (n  X  v)  --  —  f  [o5,k,v] 

c  9v 

df  (v) 

=  |  (![«»,  k]  +  v  X  B[05,k])  •  — V-  +  (B.3) 

Sv 


where  f  (0,k,v]  is  the  initial  value  of  f  (t,k,v].  Equation  (B.3) 
can  be  solved  by  the  method  of  characteristics.  By  introducing  a 
parameter  t  so  that  v  =  v(t),  the  characteristics  of  Eq.  (B.3) 


become 
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■*  t* 

V  =  as  (n  '<  T) 
c'  ' 


(a.*) 


and 


+  i(<D  -  =  g[v( 


(H  O 
\-m  •  *  r 


where  the  dot  represents  the  derivative  with  respect  to  "-he  perimeter 
t  and 


g£v]  =  |  (l[a>,k]  +  »x  B[(n,k])  •  +  ^(C,k,vj  _  (B.6) 

'  '  dv 

Note  that  t  is  not  the  "time"  which  has  been  taken  out  of  the  probloa 
by  Laplace  ecalysis.  Equation  (B,U)  is  solved  by  introducing  the  gyro- 
tensor  defined  by 

D(tp)  =  n,a  +  (I  -  n,n)c.s  9  *  tX  I  sir.  qp  (8*7} 

or,  in  tensor  rotation,  by 

Dt>>  ■  Vj  *  <‘o  *  V.)o°*  *  -  *1»  \*u’  ■ 

Then 

v(t)  =*D(-  s  o^t’)  •  v(t’)  (2.9; 

satisfies  Eq.  (B.U),  as  can  be  seen  by  inspection.  The  howogcneou*  part 
of  Eq.  (B.5)  gives 
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*•  W  ckMK*  tke  lBUfntlM  variable  to  a  =  ±  a?,  where  tke  sic?  is 
fixed  nek  that  a  >  0  for  9  >  0,  aad  ofctaia 


Ike  fvKUoa  l(£j  U  fcflM*  by  ike  i*t*fnl  of  1|.  (D2*j  «kn 
*■  C  >  3.  I*  Uw  looor  £  pin*  Z(')  and  lt«  *tt  derivative  |a(J) 
are  given  bjr  their  tae'.jtlc  coattaelloa 

X(M>(C)  -  Ix*(c*)3*  ♦  21  yr (-!)■  hB{C)«P(-C)  (d.36 

where  h  ({)  are  the  Bernite  polynoaivl*^  satisfying  the  recurrence 
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The  function  Z(£)  satisfies  the  following  differential  equation 
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To  obtain  the  relation  equivalent  to  Eq.  (0.42)  for  Im  Q  <  0,  *e  use 
Kq.  (D.36).  Continued  fraction  expansions  are  useful  in  the  intermediate 
range  of  |cj,  i.e.,  2  <  |c|  <  8,  where  neither  the  power  series 

expansion  nor  the  asymptotic  expansion  are  satisfactory. 
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Tot  «mc  Mar  tk*  ion  eye  lot  roc  and  Ion  plasat  frequencies,  the 

see*  naabers  k  are  very  snail  as  compared  to  ]k  |,  making  the  wave 
*  A 

propagation  almost  perpendicular  to  the  colunn.  In  a  certain  range  of 

anall  k  ,  |k  |  becones  so  large  that  a  prohibitive  number  of  terns 
t  A 

has  to  be  taken  la  Che  series  representation  of  the  ions  susceptibility. 


Ve  see  in  Appendix  D  that  only  one  integration  lc  needed  to  evaluate 
the  susceptibility  tensor  eleawnts,  i.e.. 
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This  integral  represents  a  Laplace  transform  of  a  product  of  two 
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11.  abstract  This  report  deals  with  the  study  of  electrostatic  waves  in  hounded  hot 
plssmss  and,  except  for  a  few  dc  measurements,  is  theoretical. 

The  study  is  broadly  divided  into  two  parts.  The  first  section  investigates 
waves  in  a  diode  system  with  no  dc  magnetic  field  applied.  A  dc  analysis  of  the 
equilibrium  plasms  produced  within  a  thermionic  diode  is  undertaken  to  provide  a 
basis  for  the  rf  analysis,  and  expected  plasms  densities  in  the  diode  are  computed 
for  a  variety  of  practical  situations.  Electrostatic  wave  resonances  in  the  diode 
are  predicted  by  using  the  hydrodynamic  model  for  the  nonunifom  plaaaa.  The  irnpe- 
dance  of  a  uniform  plssms  diode  is  obtained  by  using  a  kinetic  model  of  the  plaaaca. 
This  model  enables  U3  to  take  into  account  the  end  plate  electron  absorption  loss, 
a  process  similar  to  that  which  causes  end  plate  diffusion  in  the  Q-machine.  The 
sbsorption  loss  is  found  to  have  s  large  effect  on  the  impedance  of  the  diode. 

The  second  psrt  of  the  report  deals  with  the  study  of  guided  waves  along  a 
cylindrical  column  of  Msxwellisn  plssma  in  a  magnetic  field.  A  dc  study  of  the 
plasma  column  is  first  conducted.  Theoretical  density  and  current  profiles  are  ob¬ 
tained  snd  are  compsred  with  the  measured  results.  Since  an  rf  analysis  using  a 
self-consistent  dc  solution  is  quite  involved,  the  plasma  column  is  approximated  by 
a  uniform  plssms  with  s sharp  boundary  snd  no  drift.  To  obtain  the  rf  fields  in 
such  s  column,  a  plane-wave  solution  for  sn  infinite  Maxwellian  plasms  in  an  applied 
magnetic  field  is  first  obtained.  By  superimposing  infinitely  many  plane  waves, 
the  fields  within  the  column  are  constructed,  snd  by  matching  the  appropriate 
fields  at  the  boundary,  s  dispersion  relation  is  derived. 
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The  solutions  to  the  dispersion  relation  reveal  the  existence  of  a 
new  type  of  unstable  wave.  When  only  the  electrons  are  asstaaed  to  respond 
to  the  electric  field,  the  surface  waves  which  propagate  when  o>ce  <  cDp^ 
are  found  to  be  unstable.  When  the  ion  notion  is  included,  additional 
unstable  surface  waves  are  obtained  whenever  aid  <wpi.  A  study  of  the 
instability  shows  that  it  is  due  to  finite  Lsrnor  radii  effects  and  is 
driven  by  the  transverse  energy  of  the  particles. 

While  the  "electron"  surface  wave  instability  is  relatively  week 
and  can  easily  be  stabilised  Ly  asking  coce  >  <upe>  the  "ion"  surface 

wave  instability  is  found  to  be  very  strong,  and  the  stcbilizlng  condition  1 

of  wci  >  tDpj  is  not  easily  attained  in  high  density  plas>ias  such  ss 
those  in  fusion  machines.  Thus,  the  unstsble  ion  surfsce  waves  may  have 
a  serious  effect  on  the  containment  of  the  fusion  plasma. 
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